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ABSTRACT 

We report VLA H i Zeeman observations toward the W49A star-forming region and 
the SNR W49B. Line of sight magnetic fields (Bios) of 60 to 300 /iG at 25" resolution 
were detected toward W49A at velocities of ~ 4 km s~^ and ~ 7 km s~^. The Bi^s 
values measured toward W49A show a significant increase in field strength with higher 
resolution especially for the ~ 4 km s~^ Hi component. The Hi gas in the velocity 
range —5 to 25 km s^^ toward W49A shows good agreement both kinematically 
and spatially with molecular emission intrinsically associated with W49A. Based on 
comparisons with molecular data toward W49A, we suggest that the ~ 4 km s~^ 
H I component is directly associated with the northern part of the H ll region ring, 
while the ~ 7 km s~^ Hi component seems to originate in a lower density halo 
surrounding W49A. We estimate that the W49A North core is significantly subvirial 
(2T/I W |~ 0.2), and that the total kinetic -|- magnetic energies amount to less than 
1/3 of the total W49A North gravitational energy. These magnetic field results suggest 
that W49A North is unstable to overall gravitational collapse in agreement with 
evidence that the halo is collapsing onto the W49A North ring of Hii regions. 

The majority of the Hi column density toward W49B comes from Sagittarius Arm 
clouds along the line of sight at ~ 40 km s^^ and ~ 60 km s~^. No significant magnetic 
fields were detected toward W49B. Comparison of the spectral distribution of Hi gas 
toward W49A and W49B suggests that evidence placing W49B 3 kpc closer to the sun 
(i.e. at 8 kpc) than W49A is quite weak. Although we cannot place W49B at the same 
distance as W49A, we find the morphology of a ~ 5 km s~^ Hi component toward the 
southern edge of W49B suggestive of an interaction. 

Subject headings: ISM:clouds — ISM:individual (W49A, W49B) — ISM:magnetic 
fields — radio lines:ISM 

1. INTRODUCTION 

The W49 complex is composed of one of the most luminous star forming regions in our galaxy 
(W49A) and a relatively young supernova remnant (W49B) 12' to the east (see Fig. |l]). At high 
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resolution, W49A is resolved into numerous compact Hii regions (c.f. Dreher et al. 1984; De 
Pree et al. 1997). The Hii regions to the north (W49A North) form a "ring" of small (~ 5") 
ultra-compact H ii regions which appear to have formed during a period of triggered star formation 
(c.f. Dickel &: Goss 1990; also see Fig. ^ for a magnified view of W49 North). The Hii region with 
the highest emission measure G0, is associated with a molecular outflow and some of the strongest 
H2O masers in the galaxy (Scoville et al. 1986; Walker, Matsakis &: Garcia-Barreto 1982). To the 
SE of W49A North is W49A South which has a rim-brightened cometary shape and may have 
an edge-on interface with dense molecular gas to the west (Dickel & Goss 1990; also see Fig. |2|). 
Gwinn, Moran, &: Reid (1992) recently estimated the distance to W49A to be 11.4 kpc from H2O 
maser proper motions. The giant molecular cloud associated with the W49A star forming region 
(as traced by ^^CO) has a total mass of several 10^ Mq (Mufson & Liszt 1977). 

Efforts to explain the complex morphology and kinematics of W49A have been abundant, 
and although progress has been made, thus far no single model can reproduce its observed 
properties. Low excitation molecular emission lines (i.e. ^^C0(1 — 0)) toward W49A are resolved 
into two components with velocities of ~ 4 km and 12 km s^^ (Mufson & Liszt 1977). Such 
observations led to the suggestion that there are two distinct clouds in the vicinity of W49A, 
and that their collision is responsible for triggering the observed star-formation (e.g. Mufson & 
Liszt 1977; Miyawaki, Hayashi, & Hasegawa 1986; Serabyn, Glisten, & Schulz 1993). However, 
observations of high excitation (i.e. ^^C0(7 — 6)) and low abundance (i.e. ^^CO) lines have 
revealed the presence of colder, lower density gas at ~ 7 km s^^ which appears as self-absorption 
in more abundant species (Jaffe, Harris, & Genzel 1987; Phillips et al. 1981). Moreover, inverse 
P Cygni profiles with absorption at ~ 12 km s~^ have been observed in high resolution (~ 5") 
observations of HCO"*", NH3, and CS toward region G which indicate the presence of infall (Welch 
et al. 1987; Jackson & Kraemer 1994; Dickel et al. 1999). These facts have been taken to imply 
that the infalling gas at ~ 12 km s~^ originates from a "halo" component at ~ 7 km s~^ and is 
falling onto the the ring of H 11 regions in W49A North. Attempts to measure rotation in this gas 
or of the H 11 regions themselves have been largely unsuccessful or unreliable (Welch et al. 1987; 
Jackson &; Kraemer 1994; De Pree, Mehringer, & Goss 1997; Dickel et al. 1999). In addition to 
infall, there is also some evidence for the existence of more than one cloud or a velocity gradient 
in W49A (c.f. Dickel et al. 1999). 

W49B is estimated to lie at a distance of ~ 8 kpc (e.g. Moffett & Reynolds 1994; also see 
§4.2) and has one of the highest 1 GHz SNR surface brightnesses in the Galaxy. It has a composite 
morphology with edge brightened radio synchrotron emission but centrally concentrated thermal 
X-ray emission (c.f Moffett & Reynolds 1994; Smith et al. 1985). The spectral index of W49B 
is estimated to be a ~ —0.5 {S oc i^") and the X-ray emission is thought to originate from a 
reverse shock in the remnant's interior (Green 1988; Kassim 1989; Moffett &: Reynolds 1994). 
The presence of such a reverse shock likely indicates that this remnant is young, with an age of 
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~ 6000 years (Smith et al. 1985). MofFett & Reynolds (1994; and references therein) find that the 
synchrotron emission of W49B has an unusuahy low fractional linear polarization, despite efforts 
to detect it at high frequencies and resolutions. These authors note that the low degree of linear 
polarization observed toward W49B may be indicative of tangling in the magnetic field lines or 
Faraday depolarization within the SNR. 

Diffuse atomic and molecular gas at LSR velocities of ~ 40 and 60 km s~^ which does not 
appear to be physically associated with W49 has also been observed toward W49A and W49B. 
Some of the species observed in absorption include Hi, OH, H2CO, HCO+ and CS (Lockhart 
& Goss 1978, Sieging et al. 1982; Nyman 1984; Crutcher, Kazes, & Troland 1987; Miyawaki, 
Hasegawa, & Hayashi 1988; Greaves &: Williams 1994). The excitation requirements of these 
species suggest molecular hydrogen gas densities in the 10^ to 10^ cm~^ range. The ~ 40 and 60 
km s~^ components are thought to arise from Sagittarius spiral arm clouds. Since this line of sight 
passes through the Sagittarius arm twice, it is difficult to determine whether these clouds lie at 
their near or far distances (e.g. Dame et al. 1986; Jacq, Despois, & Baudry 1988; also see §4.2). 

Given the wide range of physical conditions that can be sampled by a single observation 
toward W49 - the complex kinematics of W49A, the SNR W49B, and diffuse gas at ~ 40 and 60 
km s~^ coupled with the limited resolution of previous W49 Hi observations (2'; Lockhart & Goss 
1978) we have performed VLA^ H i Zeeman observations toward the W49 complex. In particular, 
we hope to shed some light on the physical processes governing the star formation process(es) in 
W49A by imaging the magnetic field strengths and morphology in this intriguing source. Previous 
VLA H I Zeeman magnetic field observations toward star forming regions like W3 (Roberts et al. 
1993) and M17 (Brogan et al. 1999) have revealed line of sight field strengths of several hundred 
fiG in gas intrinsic to these sources, as well as, complex field morphologies. 

The results of our VLA H i Zeeman absorption observations toward the W49 complex are 
presented below and are organized in the following order: details of the data reduction and 
observing parameters are given in §2; discussion of the 21 cm W49A and W49B continuum 
morphology can be found in §3.1, while the properties of the Hi absorption toward W49A and 
W49B are described in §3.2; our Hi Zeeman magnetic field results are presented in §3.3; and 
finally, the implications of these data are discussed in §4. 

2. OBSERVATIONS 

The H I Zeeman data reported here consist of 12 hours of VLA B-config. and 7 hours of VLA 
D-config. data. The key parameters for these observations are given in Table |^. We observed 
both senses of circular polarization simultaneously. Since Zeeman observations are very sensitive 
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to small variations in the bandpass, a front-end transfer switch was used to switch the sense of 
circular polarization passing through each telescope's IF system for every other scan. In addition, 
we observed each of the calibration sources at frequencies shifted by it 1.2 MHz from the Hi rest 
frequency to avoid contamination from Galactic H i emission at velocities near those of W49. 

The AIPS (Astronomical Image Processing System) package of the NRAO was used for the 
calibration, imaging, cleaning, and calculation of optical depths. The right (RCP) and left (LCP) 
circular polarization data were calibrated separately and later combined during the imaging 
process to make Stokes I=(RCP + LCP)/2 and Stokes V=(RCP — LCP)/2 data sets. Bandpass 
correction was applied only to the Stokes I data sets since bandpass effects subtract out to first 
order in the Stokes V data. These H i data were also Hanning smoothed during imaging to improve 
their signal to noise (S/N). This means that while the channel separation is 0.64 km s~^, the 
velocity resolution is only ~ 1.3 km s~^. 

The AIPS task IMAGR was used to create "cleaned" Hi line+continuum data sets at four 
different resolutions: 5", 10", 15", and 25". The 25" resolution images were also convolved 
to a resolution of 40". Note that uniform weighting of these data produces a beam of ~ 5", 
while natural weighting produces images with ^ 25" resolution. Data at different resolutions 
were needed for the wide range of analysis presented here. For example, the 5" data were used 
exclusively to create high resolution 21 cm continuum images, the 10" W49A Hi line data are 
compared to H2CO and C^^O data with similar resolutions, while 15" is the highest resolution 
for which were able to make positive magnetic field detections due to S/N constraints. The 40" 
data were used to search for magnetic field detections with the greatest sensitivity. The RMS 
noise characteristics and flux density to brightness temperature (T;,) conversion factors for each 
resolution are summarized in Table ^. 

Separate line and continuum data sets were created by estimating and removing the continuum 
emission in the image plane using IMLIN for W49A and W49B separately. In this case continuum 
subtraction in the image plane turned out to be the best option for two reasons: (1) a greater 
number of line-free channels could be identified for W49A and W49B individually; (2) given the 
strength of the Hi absorption lines toward W49 (saturated at many positions), better S/N was 
obtained for the line data when the line-|-continuum were cleaned together (i.e. where the line is 
strong very little cleaning is needed - reducing cleaning errors from channel to channel). 



3. RESULTS 

3.1. 21 cm continuum of the W49 complex 

Figure [l| shows our VLA 21 cm continuum image of W49A (Hii region complex) and W49B 
(SNR) at 25" resolution. This image reveals the main continuum features of W49A and 
W49B including the numerous individual H 11 regions present toward W49A and the "blow out" 
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morphology of SNR W49B toward the NE. Inside the 15 mJy beam~^ contour level on Fig. |T|, 
the flux density from W49A is ~ 24 Jy, while the contribution from W49B is ~ 17 Jy. The rms 
noise in this 25" image is ~ 3 mJy beam~^. Evidently, about half of the 21 cm continuum flux is 
missing from our synthesis images since single dish estimates of the flux densities of W49A and 
W49B at 21 cm are ~ 45 Jy and ~ 25 Jy, respectively (Pastchenko &: Slysh 1973; Mezger et al. 
1967). Higher resolution continuum images of W49A and W49B are discussed in detail below. 



3.1.1. W49A Continuum 

Continuum images of W49A with 5" resolution at 21 cm and 6 cm are presented in Figures |^a 
and b. The 6 cm continuum data were obtained from the VLA archive and were originally 
reported in Dickel & Goss (1990). At this resolution, the wide array of distinct Hii regions which 
make up the W49A complex become apparent (see also Dickel &: Goss 1990; De Pree et al. 1997). 
A number of these H ii regions are indicated on Figure || by their letter designations following 
the nomenclature of De Pree et al. (1997). Regions A-L make up the "ring" of Hii regions first 
described by Welch et al. (1987) (see §1) and are indicated by the largest letter symbols in Fig. || 
(due to space limitations only the regions prominent at 21 cm are labeled). The Hii regions that 
are most frequently mentioned by name in the text also have arrows pointing toward them. Also 
note that with the exception of regions W49A South, LL, S, and Q the region shown in Fig. |2| is 
commonly referred to as W49A North. 

The most noticeable difference between the 21 cm and 6 cm images presented in Fig. |2| is the 
lack of extended structure present in the 6 cm image. However, this deficit is simply the result of 
a lack of short spacing information in the 6 cm VLA B-array data. The most significant difference 
between the continuum emission at these two wavelengths is that the 21 cm continuum emission 
peaks at Hll region L, as opposed to region G in the 6 cm image. This apparent lack of agreement 
is due to the fact that most of the H ii regions on the western side of the "ring" are optically thick 
at 21 cm. Indeed, Dickel & Goss (1990) found that these western Hii regions are even somewhat 
optically thick at 6 cm from comparison with data at 2 cm (also see De Pree et al. 1997). 



3.1.2. W49B Continuum 

Figures ^a and b show 5" resolution continuum images of W49B at 21 cm. The continuum 
contours on Fig. ^a reveal the small scale synchrotron emission present in the SNR, while the 
greyscale of Fig. 3b emphasizes the low surface brightness filamentary and arc like structures 
present toward the eastern and northern portions of W49B. These filamentary features are also 
observed in the multiwavelength W49B continuum images of Moffett &: Reynolds (1994). Like 
these authors, we also note the arc-like appearance of the central and western-most features along 
with the somewhat helical or twisted appearance of the eastern-most filamentary structure (also 
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see Fig. ^ in MofFett &: Reynolds 1994). Moffett & Reynolds (1994) rule out a thermal origin for 
these filaments based on high resolution (15") spectral index maps created from 90 cm, 21 cm, 
and 6 cm data. 



A number of the Hi components toward W49A (G43. 2-0.0) and W49B(G43.3-0.2) suffer 
from severe saturation effects, as might be expected from the low galactic latitude and distance 
to these sources (~ 8 to 11 kpc). Since saturated channels have an undefined optical depth 
[r = — ln(l + Tl/Tc)] (where Tl is negative and is positive; see also Roberts et al. 1995), 
saturation can cause significant underestimates of Hi column densities. This is particularly true 
of W49A and W49B, because saturated pixels are present even at the continuum peaks for some 
velocity components. To mitigate this problem we have used 



to estimate the largest optical depth value that can be reliably measured at each pixel given the 
continuum strength (Tc) and a. For the 15" (used for W49A) and 25" (used for W49B) data, 
3(T ~ 10 mJy beam~^. This "limiting optical depth map" was used to replace values in the 
original optical depth cube wherever the calculated thi is greater than the limiting value (i.e. the 
calculated value is unrealistically large because | T/^ |~ Tc) or the calculated value is undefined due 
to complete line saturation. Note that the optical depths in the channels replaced in this manner 
are lower limits to the true values. In addition, the final optical depth cubes were also masked 
wherever the continuum brightness is less than 50 and 30 mJy beam~^ for the 15" and 25" cubes 
respectively. 

While the exact values of the optical depths derived from these "corrected" optical depth 
cubes retain a degree of uncertainty (i.e. they are lower limits), we feel that they better represent 
the distribution of the Hi gas. However, two caveats should be kept in mind while examining the 
Nhi/Ts images presented in §3.2. First, despite the obvious advantage of replacing the saturated 
optical depths with a lower limit (instead of zero), this procedure can introduce a morphological 
bias which causes the Nni/Tg images to artificially resemble the continuum if large regions of 
the source at a given velocity have been replaced and also where the continuum is weak. We 
have tried to minimize this bias by identifying where this effect is present in the sections below. 
Second, although we present estimates for the Hi spin temperature Tg wherever possible (as 
deduced from previous observations), this is essentially an unknown quantity, and need not be 
constant positionally within a single component or spectrally for components of similar velocity 
(i.e. the Hi "groups" described below). Therefore, Nni/Tg images summed over multiple velocity 
components, in particular, may be biased toward components with low Tg. 



3.2. Hi Optical Depths 
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3.2.1. Hi Optical Depths Toward W49A 

The molecular gas intrinsically associated with W49A lies in the velocity range —5 to 25 
km s~^ and is kinematically complex (see eg. Dickel et al. 1999; §1). There appear to be 5 
different heavily blended Hi components in this velocity range toward W49A, but no Gaussian 
fitting was attempted due to the complexity of the spectra. The approximate center velocities 
of these Hi components are 4, 7, 10, 14, and 18 km s~^ (see for example Fig. 5). W49A Hi 
optical depth channel images in the velocity range to 21 km s^^ with 15" resolution are shown 
in Fig. ^. For comparison of the optical depth replacement scheme, the top set of panels in Fig. 
4 show thj masked wherever the displayed channel was saturated or the calculated value exceeds 
the limiting value (Eq. 1), but without any replacement, while the bottom panels have the same 
masking but saturated and uncertain optical depths have been replaced with lower limits. In both 
sets of panels, every third channel is displayed. 

Figure 4 shows that the gas near ~ 4 km s~^ is concentrated toward the ring of Hii regions 
and toward the SW (regions S and Q). In contrast, the ~ 7 km s~^ gas is widely spread over all of 
W49A, with a peak toward region J. The ~ 10 km s^^ Hi gas is almost as widely distributed, but 
shows a thi minimum between the eastern and western sides of W49A. The ~ 14 km s~^ gas has 
a distribution similar to that at 10 km s~^, but has an additional peak towards W49 South. The 
lack of H I gas between the W49 North continuum peak and the SE regions of W49A suggests that 
gas between 10 - 16 km s~^ may arise from two spatially distinct clouds. The Hi gas near ~ 18 
km s~^ has a morphology that is similar to that at ~ 7 km s~^ (i.e. widely distributed over the 
whole source). It is also noteworthy that while these velocity components have somewhat different 
overall morphologies, all of them peak toward the western side of the W49A ring of Hii regions. 
The W49A optical depth morphologies between —5 to 25 km (Fig. show overall agreement 
with the 2' resolution Hi optical depths images presented by Lockhart &: Goss (1978). 

Since these Hi data were observed in absorption, they can in principle lie anywhere in front 
of the source along the line of sight. Therefore, comparison with molecular tracers in the velocity 
range —5 to 25 km s~^ (which are intrinsically associated with W49A) are needed to establish 
whether the Hi gas in this velocity range is also directly associated with this star-forming region. 
Figure ^ shows a 10" resolution 21 cm continuum image (dashed contours and greyscale) with solid 
integrated C^^O (2 - 1) emission contours superposed (summed from —5 to 20 km s ^; Dickel et 
al. 1999). Surrounding this image are thi, 6 cm (10 x th2Co)i and (O.lx C^^O) profiles toward 
several W49A Hil regions. The VLA 6 cm H2CO data presented here have 10" resolution and 
were orig inally discussed in Dickel & Goss (1990). Similarly, the 12" resolution IRAM C^^O (2 - 1) 
data are also presented in Dickel et al. (1999) and were obtained from the UIUC Astronomical 
Digital Image Library. 

Figure |6| shows the 10" resolution Nhi/Ts integrated from —20 to 24 km s^^ (with saturated 
and uncertain optical depths replaced with lower limits as describe in §3.2). The morphology 
of the Nhi/Ts image shown in Fig. ^ is quite similar to the distribution of C^^O (2 — 1) within 
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the field of view observed by Dickel et al. (1999; see Fig. ^) in that they both peak toward the 
western side of the ring of Hii regions. The resemblance of Nni/Tg to the continuum toward the 
low surface brightness regions of W49A is a result of the saturation replacement method. The 
thi profiles shown in Fig. |5| also have a velocity extent that is quite similar to the molecular gas 
tracers H2CO and C^^O indicating that they most likely arise from similar regions in W49A. One 
exception is the Hi gas near ~ 18 km s^^ where there appears to be little molecular gas toward 
any of the W49A Hii regions (see Fig. P). Therefore, excluding the ~ 18 km s~^ gas, estimates 
of the spin temperatures (T^) of the Hi components between —5 to 17 km s~^ can be taken from 
previous dust and molecular line data of W49A. 

One reasonable estimate for the Hi spin temperature (T^) is the average C0(2 — 1) kinetic 
temperature of 51 ± 5 K estimated by Phillips et al. (1981). The ~ 20" JCMT 450, 800, and 
1100 /im images of W49A by Buckley & Ward-Thompson (1996) show that the dust emission 
peaks slightly north of region G, with secondary peaks toward the SW (regions S and Q), the NE 
(region JJ), and W49A South. The NE dust peak is positionally coincident with the enhancement 
in Thi observed toward the eastern side of W49A between 10 and 16 km s~^ (see Fig. ^). Buckley 
&; Ward-Thompson (1996) derive dust temperatures toward W49A North of 17 K, 50 K, and 145 
K; the SE of 57 K and 210 K; and SW of 48 K and 240 K (also see Sievers et al. 1991). These 
authors suggest that the cold dust toward W49 North (17 K) originates from a halo surrounding 
W49A, while the warm dust at ~ 50 K is associated with the compact Hii regions. The hot dust 
between ~ 150 to 250 K is thought to be be a minor constituent of the total dust mass. C0(7 — 6) 
observations also indicate the presence of a warm gas component (80 ^ T ^ 150) toward W49A 
North extending ~ 2', centered slightly north of Hii region G (Jaffe et al. 1987). An upper limit 
to the Hi Ts can also be estimated from the lowest 21 cm continuum brightness for which we 
can still observe Hi absorption ~ 150 K (25"), since the background continuum source must be 
warmer than the foreground absorbing gas. It seems likely then, that the Hi spin temperature in 
the -5 to 17 km s"^ range is 20 ^ ^ 150 K. 

Figures ^ and b show Nfji/Tg images for the velocity ranges 25.5 to 47.4 km s~^ and 48 to 
82 km s^^. The Hi components in these two velocity ranges are thought to arise from Sagittarius 
Arm clouds (§1). The dominant Hi velocity components in the ~ 40 km s^^ group are at ~ 34 
and ~ 40 km s~^, while the ~ 60 km s^^ group consists primarily of three components at ~ 54, 
~ 58, and ~ 64 km s~^. Overall, Nni/Tg for the ~ 40 km s~^ and ~ 60 km s~^ Hi components 
is fairly uniform over W49A. The apparent concentration of somewhat higher Nni/Tg for the 
~ 40 km s~^ group toward the eastern regions of W49A (W49A South, MM, KK) is real, as 
is the apparent increase of the ~ 60 km s^^ Nhi/Ts toward the NW. However, these regions 
also suffered from the most saturation as indicated by the black boxes in Fig. 7, and the strong 
resemblance of these peaks to the W49A continuum morphology is due to the saturated optical 
depth replacement scheme. 

It is clear from comparison of the greyscale ranges in Figs. |^, and b, that the ~ 60 km s~^ 
group has the greater overall Nni/Tg. This is because there are more Hi components in the ~ 60 
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km Hi group (i.e. greater velocity extent) since the individual Hi lines in both groups have 
approximately equal line depths. The individual Hi components in both the ~ 40 km s~^ and 
~ 60 km s~^ velocity ranges have linewidths of AvpwHM ~ 5 km s~^. As described above, a 
reasonable upper limit for Tg can be obtained from the lowest continuum brightness for which 
we can still accurately measure thi- In the case of W49A, with 15" resolution, this limit is 50 
mJy beam~^ or 150 K. 

As mentioned in §1, gas at these velocities has been observed previously in numerous 
molecular species in absorption against W49A. A recent study of these Sagittarius spiral arm 
clouds using CS(2 - 1) and CS(3 - 2) absorption by Greaves & Wilhams (1994) with ~ 20" 
resolution show that the velocities of the CS components at ~ 40 km and ~ 60 km s~^ 
are quite similar to those observed in Hi, although the CS line widths (~ 1 to 2 km s~^) are 
about three times narrower. These authors suggest that the average densities of these clouds are 
10^ ^ n ^ 10^ cm~^ for 10 ^ T^m ^ 50 K based on the relative intensities of the CS transitions. 
Additionally, Greaves & Williams find that the ~ 40 km s^^ A^(CS) is higher than that of the 60 
km s~^ CS gas, in contrast to the higher Nni/Tg found for the ~ 60 km s~^ group in Hi. The 
~ 40 km s^^ components are also measured to have the greater column density in species like CN, 
NH3, and SO; while the ~ 60 km s~^ components have greater OH, H2CO, and HCO"*" column 
densities (Greaves & Wilhams 1994; Tieftrunk et al. 1994; Nyman & Millar 1989). Evidently 
the column densities measured for these clouds are sensitive functions of resolution, abundance, 
and/or density. 



3.2.2. Hi Optical Depths Toward W49B 

The H I profiles toward W49B reveal the presence of four distinct "groups" of components at 
~ 60, 40, ~ 20, and ~ 5 km s~^. Each group consists of several Hi components that are heavily 
blended with each other, but are clearly separated in velocity from the other groups. Figures ^a 
and b show Hi optical depth profiles toward the eastern and western 21 cm continuum peaks of 
W49B and Figures ^a, b, c, d, and e show Nhi/Ts images for the total W49B Hi velocity range 
and the ~ 60, ~ 40, ~ 20, and ~ 5 km s^^ Hi groups. The 25" resolution data were used for these 
analyses (and figures) to improve the S/N of the data. We note that the unreplaced column density 
images (see §3.2) in these velocity ranges have very similar morphologies, but the magnitude of 
Nhi/Ts at most pixels is significantly less, as expected if the majority of the replaced channels lie 
near line centers. 

The first thing that is obvious from these column density images is that the morphology of 
Nhi/Ts does not particularly follow the continuum, unlike the case for the gas associated with 
W49A (e.g. Fig. 6). Another striking feature of these images is the "spotty" nature of the 
Nhi/Ts morphology. This is particularly noticeable in the ~ 20 km Nni/Tg image (Fig. 
9d). The cause of this "spottiness" is unkown, but could be indicative of unresolved structures 
since the individual spots are about the size of the synthesized beam (25"). An upper limit on 
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their corresponding linear sizescale cannot be easily determined since the distance to each of the 
Hi groups is uncertain (ranging from ^ 0.5 to 11 kpc). The total Nni/Tg image toward W49B 
(Fig. ^a) shows an increase in this ratio toward the western and SW sides of the remnant. This SW 
region of high H I column density is coincident with the region where Lacey et al. ( 2000) observe 
free-free absorption of the W49B 74 MHz continuum emission. The kinematics, morphology, and 
degree of H i saturation for the individual groups are described below. We note that the unreplaced 
column density images in these velocity ranges have very similar morphologies, but the magnitude 
of Nhi/Ts at most pixels is significantly less, as expected if the majority of the replaced channels 
lie near line centers. 

The 60 km s~^ group is composed of at least 5 components and is kinematically complex. The 
r^-i 60 km s~^ Nhi/Ts is shown in Fig. |9|b summed from 51.3 to 80.2 km s~^. With the exception 
of an Hi component at ~ 74 km s~^ which is isolated to the NE, most of the Hi gas in the ^ 60 
km s~ group is concentrated toward the western side of the remnant. Prom comparison of Fig. |b 
with Fig. it is clear that the ~ 60 group dominates the total Nhi/Ts toward W49B. Downes & 
Wilson (1974) observed an H134a radio recombination line toward W49B at a velocity of 65ibl0 
km s~^, indicating the presence of a partially ionized medium or low surface brightness Hii region 
at this velocity toward W49B. Given the similarity of the highest column density ~ 60 km s~^ Hi 
gas to the regions where low frequency free-free absorption was observed by Lacey et al. (2000), it 
seems likely that much of the H I gas in this group is intrinsically associated with the absorbing 
medium. 

The ~ 40 km s"^ Nhi/Ts image shown in Fig. ^c is summed from 26.1 to 50.6 km s ^. The 
two main components that make up the ~ 40 km s~^ group are at ~ 41 km s~^ and ~ 34.5 km s~^. 
The component at ~ 41 km s~^ is saturated toward the southern half of the remnant, while the 
~ 34.5 km s~^ component is only saturated toward the northern part of the SNR. Despite this 
saturation, it is clear that the ~ 40 Nhi/Ts has an overall concentration toward the southern half 
of the remnant. There is also a region of enhanced Nhi/Ts on the NW edge of the remnant which 
is due to a component at ~ 45.5 km s~^ that is only present toward this NW Hi concentration 
(this component is not saturated). 

The Hi group at ~ 20 km s~^ is comprised of two components and the Nhi/Ts image shown 
in Fig. ^ is summed from 13.3 to 25.5 km s^^. The Hi component at ~ 21.0 km s^^ has a fairly 
uniform column density and is not saturated. The slight enhancement of this component to the 
SE does not appear to be real. In contrast, the other component in this group (at ~ 16.5 km s~^) 
is saturated to the north but also shows a significant column density enhancement to the north. 

Perhaps the most intriguing Hi group toward W49B lies at ~ 5 km s^^ due to the coincidence 
of this velocity to those found for the molecular gas intrinsically associated with W49A (12' to 
the west). The ~ 5 km s'^ Nhi/Ts image is shown in Pig. ^ and is summed from —4.8 to 12.6 
km s~^. The 5 km s~^ components are concentrated toward the extreme southern edge of the 
remnant. This ~ 5 km s~^ Hi enhancement is coincident with a steep gradient in continuum 
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brightness which could be indicative of an interaction between the SNR shock and the ambient 
medium. The possible association of the ~ 5 km s~^ W49B Hi group with molecular gas at similar 
velocities in W49A is unclear. W49B is believed to lie some 3 kpc closer to the sun than W49A 
(e.g. Radhakrishnan 1972), making an association between the Hii region and the SNR unlikely. 
However, this distance estimate is somewhat uncertain and will be discussed further in §4.2. 



3.3. W49A Bios detections 

Stokes I and V cubes with 15", 25", and 40" resolution were fitted for the line of sight magnetic 
field strengths (Bios) using the Zeeman equation (V= 2.8i3;oc,d//2dz/ for Hi; see Troland & Heiles 
1982) and the MIRIAD fitting program ZEEMAP. The dependence of Bios on the strength of 
Stokes V and the derivative of Stokes I makes it necessary to fit only one Hi component at a 
time. This can make it very difficult to fit heavily blended lines. The complicated nature of the 
H I spectra toward W49A restricted our Bios search to those H i components that are spectrally 
distinct (on at least one side) and therefore have unblended derivatives (on at least one side). 
After identifying which Hi velocity components are good candidates in Stokes I, the Stokes V 
cube was searched in the same velocity range for the distinctive Zeeman S shaped pattern (or half 
of it if the line is blended on one side) and then fitted over the channel range that best defined 
that component. To the ensure that the detections are real, each Bios image is masked wherever 
the 25" resolution 21 cm continuum image is less than 30 mJy beam~^, and only those detections 
with a signal to noise ratio Bios/ctbi^^ > 3 are considered significant. Since random noise spikes 
can cause spurious agreements between the derivative of Stokes I and Stokes V, we also require 
that the sizescale of a given detection region equal or exceed the synthesized beam. For these 
reasons, blanked or white regions on the following Bios images only indicate regions where one (or 
all) of these conditions are not met, and do not imply that Bios is zero or even small. 

Significant line of sight magnetic fields toward W49A were detected in two different Hi 
velocity components. One of these two Bios detection zones has a velocity of ~ 4 km s^^ and lies 
to the north of Hii region G, coincident with regions I, II, and GG, while the other at ~ 7 km s~^ 
covers much of the eastern 1/3 of W49A North and is coincident with regions O and J J (see 
Fig. |2| for Hii region identification). In this context, 'coincident' only implies that the Hii regions 
are in the same direction as the the Bios detection region. A composite image showing both the 
~ 4 km s~^ and ~ 7 km s^^ Bios with 40" resolution is shown in Figure ITOl and 25" resolution 



in Figure |ll|. Several features of these Bios images should be noted: (1) The detected Bios are in 
the range ~ 60 to ~ 300 /^G; (2) The ~ 4 km s^^ Bios is negative, while the ~ 7 km is positive 
(negative Bios indicates that the line of sight field points toward the observer); and (3) The higher 



resolution 25" Bios detections shown in Figure 11 are significantly stronger than those detected 



with 40^' resolution (Fig. 10), especially for the ~ 4 km s ^ component. This last point can be 



recognized by noticing that the highest 40" Bios contour in Fig. IC is |150| /^G while this is the 



lowest 25" Bios contour shown in Fig. 11. Indeed, although an image of the 15" Bios detections are 
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not shown (significant fits at this resolution are quite patchy due to poorer S/N), there is a small 
region of significant ~ 4 km s~^ Bi^s at this resolution and the field strengths appear to be about 
a factor of two higher yet. 

At 40" resolution, in the absence of significant line blending, Bi^s ^ 30 fiG should have been 
detectable toward the continuum peak of W49A (assumes a linewidth of ~ 6 km s~^). However, 
this lower limit for detecting Bios increases in proportion to the inverse of the continuum intensity 
(i.e. at half the 40" continuum peak, Bi^g must be greater than ~ 60 fiG for it to be detectable 
in these data). Likewise, at higher resolution the continuum intensity per beam is less (if the 
source is resolved) and therefore. Bios is harder to detect. It should also be reemphasized that line 
blending can significantly inhibit the detection of Bios independent of S/N considerations. 



Sample Bios fits are shown in Figures. 12a, b, and c for the ~ 4 km s^^ component at 40", 
25", and 15" resolution. Similarly, sample 40", 25", and 15" resolution ~ 7 km s~^ Bios fits are 
shown in Figs. |^, b, and c for a position near the center of the ~ 7 km Bios detection region, 
while Figs. |l^a, and b show 40" and 25" resolution ~ 7 km s~^ Bios fits for a position near the 
top of the ~ 7 km s"'^ Bios detection region. As described above, the need to fit one velocity 
component at a time significantly restricts the number of channels that can be fit for a particular 
component. In all three figures, the channel range used in the fit is indicated by the solid portions 
of the Stokes I and dl/di' profiles. Likewise, the position, resolution, and fitted value of Bios 
zizaBi^^ printed near the top each figure. The increase of Bios with resolution described above 
can also be seen in these sample Bios fits. Unfortunately, the continuum flux and therefore Hi 
absorption toward the northern 7 km s~^ position is too weak to detect Bios at 15" resolution. 



3.4. W49B Bios (non) detections 

Similar efforts were also made to detect Bios fields toward W49B, but were unsuccessful. We 
estimate that in the absence of significant line blending. Bios ^ 40 fiG should have been detectable 
in the 40" resolution Hi data at the W49B continuum peak (assuming a linewidth of ~ 5 km s~^). 
Therefore, our lack of W49B Bios detections is in agreement with the single dish OH 1665 MHz 
Zeeman results of Crutcher et al. (1987) who tentatively detected a Bios of ~ 20 /iG toward W49B 
at ~ 60 km s^^ in OH absorption. 



4. DISCUSSION 
4.1. Magnetic Fields in W49A 
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4-1-1- Bios at ^ 7 km s ^ 

As described in §3.3, line of sight magnetic fields of ~ 200 /iG are detected toward W49A 
North to the east of region L (see Fig. ^) at an approximate velocity of 7 km s~^ (also see 
Figs. 10 & ^). The Hi component and (-B/os detection) at 7 km has a similar velocity 



to molecular emission line peaks observed in the higher order transitions of ^'^CO, C^^O, and 
Hi3co+ which are optically thin (e.g. Mufson k Liszt 1977; Jaffe et al. 1987; Nyman 1984). 
This is different from the spectra of optically thick species (i.e. ^^CO), which show instead double 
peaked profiles at ~ 4 km s~^ and 12 km s~^ (eg. Miyawaki et al. 1986). Jaffe et al. (1987) 
suggest that the optically thin ~ 7 km s^^ molecular components originate from cool halo gas 
surrounding W49A that is self-absorbed by the more abundant species like ^^CO. Indeed, their 
radiative transfer modeling of C0(1 — 0), C0(2 — 1), and C0(7 — 6) data require that T ^ 50 
K and n^^ ^ 10^ cm~^ in order to explain the presence of ~ 7 km s~^ self-absorption in the 
C0(1 — 0) and C0(2 — 1) data along with its absence in their 30" resolution C0(7 — 6) data. 

It seems quite plausible that the Hi absorption observed at ~ 7 km s~^ could arise from such 
a halo, and it is certainly the case that the 6 — 10 km s~^ Hi optical depth channel images (Fig. |^ 
show that Hi gas is quite widespread across W49A in this velocity range, as might be expected 
from a halo component. From examination of the Stokes I data cube, it also appears likely that 
the restriction of 7 km s^^ Bios detections to the eastern parts of W49A is a consequence of 
this component being more severely blended toward other parts of the source. This point can be 



appreciated by comparing the ~ 7 km s~^ component in Figs. 13 a, b, c to those in Figs. |l2^, b, 
c. It is clear from these figures that the ~ 7 km s~^ component is not very distinct toward the 
~ 4 km s~^ detection region. 

It is difficult, however, to reconcile the presence of magnetic field strengths up to ~ 300 iiG 
in gas with n^/j ^ 10^ cm~^ if the currently widespread idea that the static and wave magnetic 
energies in molecular clouds are approximately equal (eg. Brogan et al. 1999; Crutcher 1999; 
Myers & Goodman 1988a, b). The assumption of equipartition between the nonthermal, wave 
magnetic, and static magnetic energy densities implies that 

B « OAAvNTnl^^ fiG, (2) 

where Av^t ^vpyyuM is the non-thermal linewidth in km s~^, and Up is the proton number 
density in cm~^. Assuming equipartition, B ~ 300 /uG, and Avnt ~ 6 km s~^, Eq. 2 suggests 
that n//2 ~ 1 X 10^ cm~^ - in contrast to the ^ 10'^ cm~^ suggested by Jaff'e et al. (1987). 
This disagreement worsens further if the higher resolution magnetic field results showing even 



higher field strengths (see Figs. |l^ &: 11) are considered. In addition, we are only sampling one 



component of the field, so that | B | could be even higher. 

There are a number of ways that this apparent disagreement could be resolved in favor of 
equipartition. For example, if the temperature of the halo is reduced to T ^ 25 K, the upper limit 
on the density increases to ^ 10^ cm~^ (see Jaffe et al. 1987). This result is in agreement with 
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the suggestion by Dickel & Goss (1990) that the halo has a density of ~ 5 x 10^ cm~^ as inferred 
from 2 cm and 6 cm H2CO observations. A lower temperature is also supported by the ~ 20 K 
dust temperatures estimated for the halo by Buckley & Ward-Thompson (1996) and Sievers et 
al. (1991). In addition, if the Hi gas is somehow clumped or compressed relative to most of the 
molecular gas it could originate in gas of higher density than the average molecular density. In 
fact, Dickel & Goss (1990) see evidence for H2CO clumping on the sizescale of their ~ 2" beam 
toward W49A. Conversely, equipartition may simply not exist in the halo component. 



4.1.2. Bios at A km s ^ 

Of the two Biog detections, the one at ~ 4 km s~^ appears the more closely related to the 
'ring' of Hii regions in W49A since it lies to the north of Hii region G, coincident with regions I, 
II, and GG (Figs. |, 0, & There are a number of observations that indicate the presence 

of molecular gas along the same line of sight as the ~ 4 km s~^ Bi^g detections. For example, the 
~ 4 km s^^ Bios detection region is contained within the northern portions of the C^^O (2 — 1) 
integrated intensity (see Fig. 15) and the C^^O (2 — 1) profiles shown in Fig. ^ toward this region 
(i.e. regions I and II) show strong C^^O emission at ~ 4 km s~^. Additionally, the th^co channel 
images of Martin-Pintado et al. (1985) also show a strong II2CO peak toward the ~ 4 km s~^ 
Bios region. A similar peak is also observed in the H2CO data from Dickel & Goss (1990), and the 
CS (2 — 1) data of Dickel et al. (1999). Curiously however, the velocity at which these species are 
strongest is between 8 km s~^ and 15 km s~^, not ~ 4 km (see Region II II2CO profile in 
Fig. III). However, the Hi line at 4 km s^^ does not appear noticeably stronger than the other Hi 
components in this region either. 

Another clue about the properties of the ~ 4 km s~^ gas comes from the VLA radio 
recombination line velocities observed toward the cospatial Hii regions by De Free et al. (1997). 
These authors measured H92a RRL velocities of ~ 9 km s~^ for regions I and II, while GG has 
a H92a velocity of ~ 15 km s~^. This comparison suggests that the H2CO and CS(2 — 1) gas is 
closely associated with the Hii regions, while the ~ 4 km s~^ Bios gas is blueshifted toward us. 

One candidate to explain these kinematics is an outflow with a quiescent core region. However, 
the only known outflow in W49A North emanates from Hii region G, but in the opposite sense 
(i.e. the northern lobe of the outflow is redshifted, not blueshifted; see eg. Scoville et al. 1986). 
Additionally, Dickel & Goss (1990) suggest that the northern lobe of this outflow must be well 
behind the bulk of the molecular gas based on H2CO absorption measurements. Therefore, it 
seems that the outflow emanating from region G is not associated with the ~ 4 km Hi gas. 
However, given the relatively small velocity difference between the I, II, and GG Hii regions and 
the ~ 4 km s~^ Hi gas (a separation of ~ 5 — 11 km s~^) it is entirely possible that a second 
outflow with its blue lobe pointing mostly in the plane of the sky (so that its apparent Doppler 
shifted outflow velocity is small) has gone undetected. 
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From the similarity of 2 cm and 6 cm TH2C0 profiles toward W49A (2" resolution), Dickel & 
Goss (1990) conclude that that the density of the absorbing molecular gas is fairly high with an 
average value of n/^j ~ 2 x 10^ cm~^ (excluding the halo component at 6 to 10 km s~^). Using 
Eq. 2, I B |~ 400 //G based on our 15" resolution ~ 4 km s~^ Zeeman detections (see Fig. 12), 
and a linewidth of 6 km s~^, we estimate that the average density is = 1 x 10^ cm~^. This 
comparison may indicate that Biog would rise even higher if it could be measured on smaller 
sizescales. Since we have only measured the line of sight component of the field, | B j could be 
also be higher if the field is mostly in the plane of the sky. Indeed, from Eq. 2, a proton density 
of Tip ~ 4 X 10^ cm~^ and linewidth of 6 km s~^ implies that | B |~ 1500 ^G. However, it is also 
possible that the Hi gas is not well mixed with the dense molecular gas so that using the H2CO 
density is inappropriate. In the future it may be possible to measure the density and Bios of the 
gas simultaneously with molecular Zeeman tracers like CN, SO, or CCS using the Greenbank 
telescope (GET) or ALMA. 



4.1.3. Energetics of W49A 

With the new information provided by these magnetic field detections, together with 
information provided by previous W49A observations, it is now possible to investigate the energy 
balance of W49A North. To do this, it is necessary to assume that the Hi Zeeman detections 
reported here reflect the magnetic field strengths in the bulk of the molecular gas in W49A. This 
seems quite likely given the high field strengths measured toward W49A, and the commensurately 
large gas densities implied by these data (see §4.1.1 - 4.1.2). However, this assumption is 
fundamental to the following estimates, and should be kept in mind. 

Buckley &: Ward-Thompson (1996) and Serabyn et al. (1993) independently estimate from 
dust and CS observations respectively, that there is approximately 1 x 10^ within the inner 
1' (3.3 pc) of W49A North. Using these estimates, the gravitational potential energy of the inner 
1' of W49A North (i.e. R=1.7 pc) is | W |~ 4 x lO^'' erg. By comparison, the thermal+turbulent 
kinetic energy in W49A North, assuming a linewidth of 8 km s~^, is only 2T= SMa"^ ~ 7 x 10^^ 
erg (where a is the one dimensional velocity dispersion). Therefore, since 2T/| W | ~ 0.2 < 1, 
we find that the W49A North core is subvirial in agreement with Miyawaki et al. (1986) who 
found that the virial mass of the core is an order of magnitude smaller than the observed mass. If 
true, this implies that the W49A North core is subject to collapse unless there are other means of 
support. We caution, however, that the parameter values estimated above are uncertain to a least 
factors of two and that these uncertainties could conspire to increase 2T/| W |. 

The idea that magnetic fields (both the static and wave components) can help support 
star- forming clouds against overall collapse has long been suspected (see e.g. McKee 1999). In 
the case of W49A North we have detected average line of sight magnetic field strengths (Bios) of 
~| 300 [ fiG at 25" resolution (even higher at 15" resolution) implying that a reasonable estimate 
for the total average field strength is | B |~ 600 fiG. The static magnetic energy resulting from a 
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field of this magnitude is A4s={t>/3)B'^R^ r-^ 5 x 10^^ erg (where b depends on the field geometry 
and is ~ 0.3; see McKee et al. 1999). Moreover, since the wave component of the magnetic field is 
thought to be in equipartition with the turbulent kinetic energy, there is an additional contribution 
to the magnetic energy of Ad^j'^T^ 3.5 x 10^^ erg (see e.g. Brogan et al. 1999; Zweibel Sz McKee 
1995). From these estimates, it appears that the the static and wave components of the total 
magnetic energy are not equal, in contrast to the results for a number of other regions including 
the M17SW star-forming core (Brogan et al. 1999; see also Crutcher 1999). Additionally, these 
estimates show that the total magnetic energy is only ~ 0.1| W | and that the magnetic and 
kinetic energy together only add up to about a quarter of the gravitational energy. Even if the 
total static field strength is as high as | B |~ 1500 fiG as suggested in §4.1 based on estimates for 
the density of the W49A core, A^s~ 3 x 10^^ erg and the magnetic plus kinetic energies still only 
add up to ~ 0.34| W |. 

While it must be emphasized that the parameters used in the energy calculations above are 
fairly uncertain, the gravitational energy in the W49A North core appears to be ~ 3 to 4 times 
larger than the combined magnetic plus turbulent kinetic energies. This imbalance suggests a 
natural explanation for both the copious amount of star-formation on-going in W49A and the 
apparent overall collapse of the W49A North halo onto the ring of H li regions suggested by Dickel 
et al. (1999). Indeed, ignoring the external pressure, we estimate that the static component of the 
magnetic field would need to be on the order of 5 mG to obtain | W | = 2T + A4w + A4s (using 
M ^~7'). While the increase in BiQg with resolution observed in these data suggest that tangling 



or small scale structure in the field could reduce the measured values of Bios (see e.g. Figs. IC, 11 



12, 13), and we are measuring only one component of the field, it seems unlikely that these effects 



would amount to the necessary factor of 10 difference between the highest observed Bios and | B 



4.1.4- Magnetic Field Morphology 

As shown in §3.3, the ~ 4 km s~^ and ~ 7 km s~^ Bigg detections point in opposite directions 
(the ~ 4 km s~^ Bigg points toward us), and lie toward different regions of W49A (see Figs. 10 
& 11). It is impossible to determine the significance and implications of this morphology without 
knowing the plane-of-sky magnetic field direction and magnitude, and the effect of line blending 
on detecting i?/os- However, two facts are clear from the Hi Bios data presented here: (1) while 
B 

los was only detected at ~ 7 km s ^ toward the eastern regions of W49A North, this is also the 
region where the ~ 7 km s~^ Hi gas is most distinct as a separate velocity component (see also 
§4.1.1); (2) the ~ 4 km s~^ component is quite distinct in regions where a significant Bios at this 
velocity was NOT detected (compare for example Figs. 12 & 13). From these facts, we suggest 
that the confinement of significant ~ 7 km s~^ Bios to the eastern side of W49A may simply be 
the result of line blending further west, and that the strong ~ 4 km s^^ Bios fields are restricted 
to the NW region of W49A. 

In the absence of more data, particularly about the ~ 7 km s~^ Bios toward the same regions 
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where the ~ 4 km s Bios is detected, there are many different scenarios that can explain the 
opposite signs of the ^ 4 km s~^ and ~ 7 km s~^ Bi^g in different parts of the source. For example, 
the negative direction of the ~ 4 km s~^ Bi^g could be the result of an outflow like the one 
proposed in §4.1.2 to explain the kinematics of the ~ 4 km s~^ component (which is blueshifted 
with respect to the nearby Hii regions). Alternantively, the infall of the halo gas reported by 
Dickel et al. (1999; also see §1) could also potentially play a role in the apparent field reversal. 
Further progress in understanding the magnetic field morphology of this complex source will have 
to await future linear polarization and molecular Zeeman observations. 



4.2. Comparison of H i gas toward W49A and W49B 



Figure 16 compares the average 25" resolution Hi optical depth toward W49A with that 
observed toward W49B. Before calculating the average, the optical depth cubes were masked 
wherever the 21 cm continuum flux is less than 30 mJy beam~^. The overall shapes of the average 
Tffj profiles toward W49A and W49B shown in Fig. |l^ (upper panel) are very similar to those 
presented in Lockhart & Goss (1978) using the Owens Valley interferometer with 2' resolution. 
However, the magnitudes of the average optical depth profiles shown here are smaller by about a 
factor of 1.5. This discrepancy is caused by our replacement of saturated and uncertain optical 
depths with conservative lower limits. This procedure underestimates thi, especially toward the 
weaker parts of the continuum (see §3.2), causing an overall lowering of the average optical depth. 
Since this method was applied equally to W49A and W49B it should not change the relative 
differences between them. Indeed, the similarity of the average 'W4:9A{thi)—W49B{thi) profile 
presented in Fig. |l^ (lower panel) to the one shown in Lockhart & Goss (1978) suggests this is the 
case. 

Although the distance to W49A has been accurately determined to be 11.4 kpc from H2O 
maser proper motion studies (Gwinn et al. 1992), the distance to W49B remains uncertain. Given 
that there is strong Hi absorption up to velocities of ~ 75 km s~^ for both W49A (G43.2— 0.0) and 
W49B (G43.3— 0.2), it is clear that W49B must he beyond the tangent point at ~ 6.2 kpc (using 
galactic center distance of 8.5 kpc). In addition, since the Hi absorption toward W49B does not 
extend to negative velocities, the SNR must also lie within the Solar circle at ~ 12.4 kpc. In the 
past, comparison of the average thi profiles W49A and W49B have also been used to estimate 
the distance to W49B. The major difference between the average Hi optical depths toward W49A 
and W49B is the absence of Hi absorption at velocities between ~ 7 to 14 km s~^ and ~ 50 to 
55 km s~^ toward W49B (see also Lockhart Sz Goss 1978; Radhakrishnan et al. 1972). Indeed, 
the greater overall optical depth toward W49A from ~ 7 to 14 km s^^ and in the ~ 60 km s~^ 
components has been used in the past to suggest that W49A lies farther away than W49B by ~ 3 
kpc (Kazes &: Rieu 1970; Radhakrishnan et al. 1972). 

We do not find clear evidence in the high resolution Hi data present here that this is the 
case. For example, the Nni/Tg images in all velocity ranges toward both W49A and W49B 



- 18 - 



show noticeable changes on size scales of ~ 1' (see Figs. ^ 0, and |^. Therefore, given the 12' 
separation between them, it is not surprising that there are significant variations in the average 
optical depth profiles toward W49A and W49B regardless of their relative distances. Indeed, 
although the average thi from 50 to 58 km s~^ is stronger toward W49A, there appears to be 
more Hi gas in the ~ 58 to 65 km range toward W49B (Fig. 16). Also, while there is a 
complete absence of Hi absorption at velocities between ~ 7 to 14 km s~^ toward W49B, there are 
a number of scenarios involving differences in kinematics, distribution, and temperature that could 
explain this behavior without assuming that W49B is closer. In any case, conversion of the excess 
W49A ~ 7 to 14 km s~^ Nni/Tg to the distance between W49B and W49A strongly depends 
on the assumed density and temperature of the Hi gas. For example, if the 7 to 14 km s~^ 
Hi absorption arises from the W49A halo component as suggested in §4.1, (n/^j ~ 10^ cm~^ and 
Ts ~ 20 K), then the necessary line of sight separation between W49A and W49B reduces to ~ 1 
pc compared to the 3 kpc obtained by Kazes & Rieu 1970 (also see Radhakrishnan et al. 1972). 
Summarizing, comparison of 25" resolution H i absorption toward W49A and W49B provides no 
clear evidence that W49B lies closer than W49A (i.e. at 8 kpc). 



5. SUMMARY AND CONCLUSIONS 

The Hi gas in the velocity range (—5 to 25 km s~^) shows good agreement both kinematically 
and spatially with molecular emission intrinsically associated with W49A (see Figs. ^ and |5|). 
Therefore, the Hi gas in this velocity range is likely to originate near the W49A star- forming 
complex. 

Significant line of sight magnetic fields toward W49A were detected in two different Hi 
velocity components. Some of the properties of these detections are: (1) Bios of 60 to 300 /iG were 
detected toward W49A at velocities of ~ 4 km s~^ and ^ 7 km s^^; (2) The ~ 4 km s^"^ Bios is 
negative, while the ~ 7 km s^^ is positive (negative Bios indicates that the line of sight field points 
toward the observer); and (3) The Biog values measured toward W49A show a significant increase 
in field strength with higher resolution especially for the ~ 4 km s~^ Hi component (see Figs. 10, 



11, |T2|, 13, and |T^). Based on comparisons of the ~ 4 km s^^ and ^ 7 km Bios detection 
regions and velocities with molecular data toward W49A, it seems likely that the ~ 4 km s^-*^ 
Hi component is directly associated with the northern part of the Hii region ring, while the 
~ 7 km s^^ H I component seems to originate in a lower density halo surrounding W49A. 

From previous molecular line and dust observations of W49A, we estimate that the W49A 
North core (~ 1' in extent) is significantly subvirial, with 2T/[ W 0.2. In addition, it appears 
that the total magnetic energy (assuming | B |~ 600 /iG and T) only contributes 0.1| W |. 

Adding the different energy contributions, we estimate that the total kinetic + magnetic energies 
only amount to less than about 1/3 of the total gravitational energy. Indeed, we find that 
I B |~ 5 mG is needed for there to be overall equilibrium in W49A North (ignoring external 
pressure). Therefore, while it must be emphasized that the parameters used in these estimates are 
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fairly uncertain, these new magnetic field results suggest that W49A North is unstable to overall 
gravitational collapse as implied by evidence that the halo is collapsing onto the W49A North ring 
of Hii regions (e.g. Dickel et al. 1999). 

There are four main Hi velocity groups detected toward W49B with velocities of ^ 60, ~ 40, 
~ 20, and ~ 5 km s~^. The majority of the Hi column density toward W49B comes from the ~ 60 
and ~ 40 km s~^ Hi components which arise from Sagittarius Arm clouds along the line of sight. 
No significant Bi^s fields were detected toward W49B. Comparison of the spectral distribution of 
Hi gas toward W49A and W49B (Fig. ^) suggests that there is no clear evidence that W49B is 
3 kpc closer to the sun than W49A (known to be at 11.4 kpc from water maser proper motions; 
Gwinn et al. 1992). Although we cannot place W49B at the same distance as W49A, we find 
the morphology of the ~ 5 km s~^ H I gas toward the southern edge of W49B suggestive of an 
interaction. Indeed, gas at similar velocities is known to be intrinsically associated with W49A. 
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1Q^OS"^54^ 36^ 24^ 12^ OO^ 07'^48^ 



Right Ascension (1950) 

Fig. 1. — Greyscale 21 cm continuum image of the W49 complex with the Hll region W49A to 
the west and the SNR W49B to the east. Continuum contours at -15, 15, 50, 100, 200, 300, 450, 
600, 900, 1200, and 1800 mJy beam'^ are also shown. The resolution is 26'.'5 x 23'.'9 (P. A. -50.8°). 
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Fig. 2. — W49 A continuum images with 5" resolution, (a) Greyscale 21 cm continuum image 
of W49 A witli 5" resolution. Continuum contours are shown at 10, 30, 50, 70, 90, 110, and 130 
mJy beam~^. (b) Greyscale 6 cm continuum image of W49 A with 5" resolution. Continuum 
contours are shown at 50, 100, 200, 300, 500, and 900 mJy beam~^ (original VLA data from Dickel 
Sz Goss 1990). The individual Hll regions that are bright at 21 cm are indicated using the naming 
convention of De Free et al. (1997) and are repeated on the 6 cm image for reference. The largest 
letter designations indicate the Hli regions which make up the "ring". Notice that at 21 cm the 
continuum peaks toward region L, while at 6 cm the continuum peak is toward Hll region G. 
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Fig. 3. — Greyscale 21 cm continuum images of the SNR W49 B with 5" resolution, (a) Continuum 
contour levels are 5, 10. 15, 20, 25, 30, 35, 40 mJy beam~^. (b) The greyscale intensity has been 
adjusted to emphasize the continuum filaments visible in the northern and NE parts of the SNR 
(see also Moffett & Reynolds 1994). 
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Fig. 4. — Greyscale Hi optical depth channel images toward W49A in the velocity range —0.27 
to 21 km s"^ (H I gas in this channel range appears to be directly associated with W49A). The 
top set of panels show the optical depths masked as described in §3.2 but without any replacement 
of saturated optical depths, while the bottom panels have the same masking but saturated optical 
depths have been replaced with lower limits. These H I data have 15" resolution and only every 3rd 
channel is shown. The contours show the 15" resolution 21 cm continuum at 50, 150, 450, and 750 
mJy beam~^. 
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Fig. 5. — Image of the W49A 21 cm contimmm [greyscale and dashed contours) superposed 
with C^^O (2 - 1) integrated emission contours summed from —5 to 20 km s ^ {solid contours). 
Surrounding this image are thi {thick black), 10 x TH2C0 {dotted) and 0.1 x C^^O emission {thin 
black) profiles toward several of the W49A Hll regions. The {dashed) 21 cm contours are at 30, 
130, 230, 330, and 430 mJy beam~^. The resolution of the 21 cm continuum, thi and TH2CO data 
is 10", while the C^^O data have a resolution of 12". The VLA H2CO 6 cm data presented here 




Fig. 6. — Greyscale image of the lower limit of Nhi/Ts toward W49A integrated from —20 to 24 
km with 10" resolution and white contours at (0.6, 0.8, 1.0, 1.2, 1.4, and 1.6) xlO^° cm~^ K~^. 
21 cm continuum contours (dashed) at 30, 130, 230, 330, and 430 mJy beam~^ with 10" resolution 
are also shown. 
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Fig. 7. — Greyscale images of the lower limit of Nni/Tg toward W49A for (a) the ~ 40km s^^ 
and (b) the ~ 60 km Hi components toward W49A with 15" resolution. The ~ 40 km s~^ 
Nhi/Ts is summed from 25.5 to 47.4 km s^^ and the ^ 60 km s^^ Nui/Tg is summed from 48 to 
82 km s^^. The white Nni/Tg contours are located at (0.55 and 0.65) xlO^'^ cm~^ on (a) and 
(0.8, 0.9, and 1.0) xlO^° cm~^ on (b). The black boxes superposed on (a) and (b) indicate 
regions where a significant fraction of the integrated channels were saturated and therefore replaced 
with lower limits (see §3.2). 21 cm continuum contours {dashed) with 15" resolution are also shown 
at 50, 150, 450, and 750 mJy beam~^ 
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Fig. 8. — Hi optical depth profiles toward (a) the eastern W49B 21 cm continuum peak and (b) 
the western W49B 21 cm continuum peak. These data have 25" resolution and the separation 
between the two peaks is ~ 3'. The channel near ~ 40 km in (b) is saturated and was replaced 
with a lower limit as described in S3. 2. 
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Fig. 9. — Greyscale integrated images of the lower limit of Nni/Tg x 10^'^ cm^^ for the (a) 
total Hi velocity range, (b) ~ 60 km Hi group, (c) ~ 40 km Hi group (d) ~ 20 km 
Hi group, and (e) ~ 5 km Hi group toward W49B (25") . The black contours show the 21 cm 
continuum at 0.02, 0.05, 0.1, 0.2, 0.3, 0.45, and 0.6 mJy beam"^ with 25" resolution. 
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Fig. 10. — Composite greyscale image of Bios for both the 4 km s~ and 7 km s~ Hi components 
toward W49A with 40" resolution. The black contom's show the 21 cm continuum at 40, 140, 
240, 340, and 440 mJy beam~^ with 10" resolution (same as Fig. ^). The white contours indicate 
Bios at —100 and —150 fiG for the 4 km s~^ component and +100 and +150 fiG for the 7 km s~^ 
component. 
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Fig. 11. — Composite greyscale image of Bi^g for both the 4 km s~ and 7 km s~ Hi components 
toward W49A with 25" resolution. The black contours show the 21 cm continuum at 40, 140, 
240, 340, and 440 mJy beam~^ with 10" resolution (same as Fig. ^). The white contours indicate 
Bios ^-t —150 and —250 iJ,G for the 4 km s~^ component and +150 fiG for the 7 km s^^ component. 
The increase in Biog with resolution is particularly noticeable in the 4 km s^^ component when this 
image is compared to Bios in Fig. 
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Fig. 12.— Bfo, fits toward W49A at 4 km s"^ (position is 19^07™51.0", +09°01'35") with (a) 40", 
(b) 25", and (c) 15" resolution. The upper panels show the VLA Stokes I profiles (solid histogram), 
and the bottom panels show the VLA Stokes V profiles (solid histogram) with the fitted derivative 
of Stokes I shown as smooth dotted curves. The value of Bi^s fit for each position and its calculated 
error are given at the top of each plot. The solid portion of the Stokes I histogram [upper 
shows the velocity range used in the fit. 
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Fig. 13.— Bios fits toward W49A at 7 km s"^ for position 19'^07"53% +09°01'06"with (a) 40", (b) 
25", and (c) 15" resolution. The upper panels show the VLA Stokes I profiles {solid histogram), 
and the bottom panels show the VLA Stokes V profiles {solid histogram) with the fitted derivative 
of Stokes I shown as smooth dotted curves. The value of Bios fit for each position and its calculated 
error are given at the top of each plot. The solid portion of the Stokes I histogram {upper panel) 
shows the velocity range used in the fit. 
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Fig. 14.— Bi^, fits toward W49A at 7 km s"^ for position 19'^07™57^ +09°0l'49"with (a) 40", 
and (b) 25" resolution. The upper panels show the VLA Stokes I profiles {solid histogram), and 
the bottom panels show the VLA Stokes V profiles {solid histogram) with the fitted derivative of 
Stokes I shown as smooth dotted curves. The value of Bios fit for each position and its calculated 
error are given at the top of each plot. The solid portion of the Stokes I histogram {upper panel) 
shows the velocity range used in the fit. 
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Fig. 15. — Composite greyscale image of Bi^g for both the 4 km and 7 km Hi components 
toward W49A with 40" resolution. Solid black integrated C^^O (2 — 1) emission contours with 12" 
resolution (integrated from —5 to 20 km s~^) are also shown (Dickel et al. 1999; also see Fig. |5|). 
The black {dashed) contours show the 21 cm continuum (10" resolution) at 40, 140, 240, 340, and 
440 mJy beam~^ (same as Fig. 10). 
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Fig. 16. — Upper panel: Comparison of average H I optical depth profiles toward W49A (solid) 
and W49B (dotted). Lower panel: Difference between the average Hi optical depth profiles shown 
in the upper panel [W49A(rij7)-W49B(ri//)]. 
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Table 1. Parameters of VLA Hi Zeeman Observations Toward W49 



Parameter 



Value 



Frequency 

Observing date for D array data 
Observing date for B array data 
Total Observing time in D array 
Total Observing time in B array 
Primary beam HPBW 
Phase and pointing center(B1950) 
Frequency channels per polarization 
Total bandwidth 
Velocity coverage 
Channel separation 
Velocity resolution " 
Angular to linear scale 



1420 MHz 
Mar 12, 1999 
Jan 2, 2000 
7 hr 
12 hr 
30' 

19''08™20.0^ +09°00'00" 
256 

781.25 kHz (165.05 km s^^) 
-52.5 to +112.5 km 
3.052 kHz (0.64 km s^^) 
1.8 km s~^ 
10" = 0.55 pc 



After Banning smoothing. 
^Assumes distance of 11.4 kpc. 
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Table 2. Parameters for W49 Images of Different Resolutions 



Resolution 


RMS (line) 


RMS (cont.) 




(") 


(mJy beam^^) 


(mJy beam~^) 


10 mJy beam^^ = Ti,(K) 


5 




8 


272 


10 


5 


5 


68 


15 


3 


4 


30 


25 ^ 


2 


3 


11 


40 


4 




4 



^Conversion factor to change 10 mJy beam into temperature in K for each 
resolution. 

''Exact resolution is 26''5 x 23'/9, RA. -50.8°. All other resolutions are exact 
with P.A. = 0.0° after convolution. 



